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ABSTRACT 

With the application to weld defects in mind, the interaction problem bet- 
ween a planar crack and a flat inclusion in an elastic solid is considered. 

The elastic inclusion is assumed to be sufficiently thin so that the thickness 
distribution of the stresses in the inclusion may be neglected. The problem 
is reduced to a system of four integral equations having Cauchy type dominant 
kernels. The stress intensity factors are calculated and tabulated for various 
crack-inclusion geometries and the inclusion to matrix modulus ratios, and for 
general homogeneous loading conditions away from the crack-inclusion region. 


1. Introduction 

In studying the strength and fracture of structural solids it is often 
necessary to take into' account, among other factors, the effect of the imper- 
fections in the material. Generally such imperfections are in the form of 
either geometric discontinuities or material inhomogeneities. For example', 
in welded joints, various shapes of voids, cracks, notches and regions of lack 
of fusion may be mentioned as examples for the former and variety of inclu- 
sions for the latter. From a viewpoint of fracture mechanics two important 
classes of imperfections are the planar flaws which may be idealized as cracks 
and relatively thin inhomogeneities which may be represented by flat inclu- 
sions. 

The correct way of modeling an inclusion would perhaps be to consider 
it as an elastic continuum fully bonded to the surrounding matrix. In this 
case, however, the crack- inclusion problems are generally difficult and only 
simple geometries and orientations can be treated analytically (see, for 


This study was supported by the U.S. Department of Transportation Office of 
University Research under the Contract DTRS 5682-C-00014 and by NASA-Langley 
under the Grant NGR- 39-007-011 . 
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example, [1], [2]). A simple feature of such crack-inclusion interaction 
problems is that generally the stress intensity factors are magnified if the 
stiffness of the inclusion is less than that of the matrix and are diminished 
if the inclusion is stiffer than the matrix. For certain types of "flat" 
inclusions a simpler way of modeling may be to represent them as either a 
membrane with no bending stiffness or a perfectly rigid plane stiffener with 
negligible thickness. In these problems one may use the basic body force 
solution as the Green's function to derive the related integral equations. 

On the other hand, since the flat inclusion with an elastic modulus smaller 
than that of the matrix would itself have a behavior similar to a crack, it 
needs to be modeled basically as a "cavity" rather than a "stiffener". 

Even though the technical literature on cracks, voids and inclusions 
which exist in the material separately is quite extensive, the problems of 
interaction between cracks and inclusions do not seem to be as widely studied. 
Such problems may be important in studying, for example, the micromechanics 
of fatigue and the fracture in welded joints. In this paper a simple model 
for flat elastic inclusions is presented and the crack-inclusion interaction 
problem is considered for various relative orientations. 

2. Integral Equations of the Problem 

The plane strain or the generalized plane stress interaction problem 
under consideration is described in Fig. 1. It is assumed that the boundaries of 
the medium are sufficiently far away from the crack-inclusion region so that 
their effect on the stress state perturbed by the crack and the inclusion 
may be neglected and the plane may be considered as being infinite. 

Referring to Fig. 1 we define the following unknown functions 


MV = 3 x^ C v i < x i > + °)~ v i ( x i »“°)II > (a<x 1 <b) , 

(i) 

h-j(x-|) = Cu 1 (x r+ 0)-u 1 (x r -0)] ^ (a<x-j<b) , 

(2) 

9 2 ( X 2 ) = "gx ^2^2 ^ * (c<X2 < d) * 

(3) 

h 2 ( x 2 ) ” ^2^2 5 "^^ ^2^2* » (c<X2<d) 

(4) 
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where u and v are, respectively, x and y components of the displacement vec- 
tor in the coordinate systems shown in the figure. It is assumed that the 
inclusion fills a flat cavity the initial thickness of which is h Q (x) which 
is "small" compared to its length 2a-j. It is also assumed that the thickness 
variation of the stresses and the strain zl in the inclusion are negligible. 
Thus, for the plane strain case, from the Hooke's law we obtain the fallowing 
stress-strain relations in the inclusion 


£ yy (x ^ = 


1 ‘ v o" 2v o 

w 


°yy (x l } 


£ xy (x l } = 2 u q 


a xy (x ^ 


(5) 


where E„ , v , y are the elastic constants of the inclusion. Now, by observing 
o o o 

that 

e yy (x-,) = [v 1 (x 1 ,+0)-v 1 (x 1 ,-0)]/h Q (x 1 ) , . (6) 


and 


2 £ xy( x l ) Cu 1 (x 1 ,+0)-u 1 (x 1 ,-0)]/h o (x) , 

E rt = 2y(l+v ) , k = 3-4v , 

o M o o' o o 


from (1), (2) and (5)-(8) we find 

, i X-| 

; K+\ jl- 

i / \ 0 0 ! 

cr t 


yy (x i } = <°-i h^T J 9 i (t)dt * 


i 


(7) 

( 8 ) 

(9) 


a 


i 

xy 



H o 

W 


h^ (t)dt 
a 


(10) 


If we let the medium to be uniformly loaded away from the crack-inclusion 
region as shown in Fig. 1, for the stress components along the x-| and axes 
we obtain 


l°o , » oo l 00 / n\ OT 

°yy (x l >°> = °yy ’ a xy (x l ’ 0) = a xy ’ 

0 y"(x 2 ,O) = a” y cos 2 e + °" x sin 2 e - 2a" y sinecose , 


(11) 

(12) 
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From the basic dislocation solution given in, for example, [3], referred 
to the coordinate system the stress state at a point (x^y-j) in the plane 

due to the displacement derivatives g-j ,h.| defined by (1) and (2) may be 
expressed as 

b 

a xx^ x l ,y l^ = | ^ G xx^ x l ’ y l st ^ g l + H xx^ x l ,y l ’O^OOldt > 04) 

a 

b 

ojj(x 1 ,y 1 ) = | [G yy (x 1 ,y 1 ,t)g 1 (t) + H yy (x ] ,t)h 1 (t)]dt , (15) 

a 
b 

CT xy( x l ,y l ) = { CG xy (x r y l ,t)g l (t) + H yy( x i’ y l’ t)h l (t)]dt ’ (16) 

a 

v/here 

G xx (x,y,t) = A(t-x)[(t-x) 2 -y 2 ] , 

G yy (x,y,t) = A(t-x) [3y 2 +(t-x) 2 ] , 

G (x,y,t) = Ay[y 2 -(t-x) 2 ] , 

Ay 

H xx (x, y ,t) = Ay[y 2 +3(t-x) 2 ] , (17) 

H yy (x,y,t) = Ay[y 2 -(t-x) 2 ] , 

H (x,y,t) = A(t-x)[(t-x) 2 -y 2 ] , 

Ay 


A(x,y,t) 



and y and k are the elastic constants of the medium (y=E/2(l+v), <=3-4v for 

plane strain and k=(3-v)/(1+v) for generalized plane stress). Similarly, 

o? 

referred to the axes x 2 ,y 2 the stress state (i,j=x,y) in the plane due 
to g 2 »h 2 may be obtained from (14)-(17) by substituting (c,d) for (a,b) and 
(x 2 ,y 2 ) for (x-j.y-,) and (g 2 ,h 2 ) for (g-, ,h-j ) . 
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The integral equations to determine the unknown functions g-j>h.|,g 2 , 
and h 2 may be obtained from the following traction boundary conditions along 


(y^=0, a<x-|<b) and (y 2 =0, c<x 2 <d): 

°yy ( x i j 0 ) + a yy( x r°) + °yy( x r 0 ^ = a yy( x l) » (a<x-,<b) , (18) 

°xy (x r 0) + a xy (x l’ 0) + a xy (x l’ 0) = CT xy (x l ) ’ ( a<x l <b > ’ < 19 > 

cr yy( X 2 ’0) + a yy (x 2’ 0) + a yy (x 2’ 0) = 0 » ( c<x 2 <d) » (20) 

a^y( x 2 ,°) + tfxy( x 2>°) + a xy^ x 2’°^ = 0 » ( c<x 2 <d ) » ( 21 ) 


where all except the coupling stresses in the second column are given by 
(9)-( 17) . The coupling stresses have the following meaning: CT yy( x i »0) is 

the normal stress on y-|=0 plane due to the displacement derivatives g 2 (x 2 ) and 
h 2 (x 2 ) and a 2 y(x 2 ,0) is the normal stress on y 2 =0 plane due to g-| , h-| , etc. 
Thus, after making the necessary stress transformations similar to (12) and 


(13), we obtain 

d 

0 yy (x r O) = j £Gyy(x r t)g 2 (t) H yy( x q ,t)h 2 (t)]dt , ( 22 ) 

c 
d 

oJy( x i,°) = E^xy( x l »t)g 2 (t) + H Jy(x 1 ,t)h 2 (t)]dt , (23) 

c 
b 

ayy( x 2 ,°) = [Gyy(x 2 ,t)g 1 (t) + Hyy(x 2 ,t)h 1 (t)]dt , ( 24 ) 

a 

b 

0 xy( x 2’°) = \ [G x]' (x 2> t>9 l (t) + ^xy( x 2 ’t)^l ( t ) 3 dt , (25) 

a 

where from 

<7yy( x r°) = <7yy(x 2 ,y 2 )cos 2 e + a^sin 2 e + a^ 2 sin2e (26) 


calculated at x 2 = x iCose, y 2 =-x-|Sin0 we have 
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( 27 ) 


Gyy( x i>t) = G yy (x 1 cos9,-x- | sin8,t)cos 2 9 + G X x^ cos9, " x l sin9 ’ t ) 
+ G xy^ x l cos0 ’ x i sin9 ) s i n2s , 

Hyy( x i »t) = H yy (x 1 cos9,-x-|Sin9,t)cos 2 9 + H xx (x-|Cose,-x 1 sine,t) 
+ H xy (x lC os9,-XiSin9,t)sin2e . 

Similar expressions for the remaining kernels in ( 23 ) - ( 25 ) are 
by using the stress transformations 

a xy( x l *°) = ^xx^ x 2* y 2^°yy^ x 2* y 2^ sin0cos0 

+ a xy (x 2 ,y 2 )(cos 2 9-sin 2 9) , (x 2 =x-|Cose, y 2 =-x-,sine) 

CT yy( x 2’°) = a yy (x l’ y l )cos20 + a xx (x l ,y l ) sin2e 

- a xy( x i »y ] ) s in2e , (x-,=x 2 cos9, y^XgSine) , 

a xy (x 2 ,°) = [a yy (x 1 ,y 1 )-a xx (x 1 ,y ] ) ]sin9cos0 

+ a xy (x 1 ,y 1 )(cos 2 e-sin 2 e), (x 1 =x 2 cos0, y 1 =x 2 sine) . 
Thus, from (14)-(25) and (29)-(31) it follows that 

G xy^ x l = C G xx^ x,y,t ^“ G yy^ x,y,t ^ s1n0cos0 

+ G xy (x,y,t)cos29, ( x = x-j cos9 , y=-x-|Sine) , 

H xy( x l st ) = t H X x( x>y » t ^" H yy^ x ’ y,t ^ sin0cose 

+ H xy (x,y,t)cos2e, (x=x 1 cose, y=-x-|Sine) , 


(28) 

obtained 

(29) 

(30) 

(31) 


(32) 

(33) 
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( 34 ) 


G yy (x 2’ t} 


Hyy( x 2’t) 


G xy( x 2 at ) 


= G yy (x,y,t)cos 2 9 + G xx (x,y 5 t)sin 2 e 

- G xy (x,y,t)sin2e , (x=x 2 cose, y=x 2 sine) 
= H yy (x,y,t)cos 2 e + H xx (x,y,t)sin 2 6 

- H (x,y,t)sin29 , (x=x ? cos9, y=x ? sin9) 

A Jr ^ 


[ G yy ( X 5 y 5 t )- G XX ( X ’ y > t )] Sin9COS9 


(35) 


+ G xy (x,y,t)cos29 , (x=x 2 cos9, y=x 2 sin ) , 


H xy (x 2 ,t) = [H yy (x,y,t)-H xx (x J y,t)]sinecose 

+ H xy (x,y ,t)cos2e , (x=x 2 cos9, y=x 2 sine) 


(36) 


(37) 


From (18)-(21 ) the integral equations of the problem may then be obtained as 
b x, d 


1 


u A 1 r 

g-| (t)dt + G(x-,)g 1 (t)dt + c Q J G yy (x ] ,t)g 2 (t)dt 




H yy( x r t )h 2 (t) dt = " c 0 a yy » (a<x-,<b) 


(38) 


1 

7T 


^-h^tjdt + j H(x ] )h 1 (t)dt + c Q j G xy (x ] ,t)g 2 (t)dt 


. 12 / 


+ c. 


d 

f „12 


H xy ( x 15 t) h 2 ( t ) dt = _ c 0 a xy ’ ( a<x i <b ) » 


(39) 


u D U 

G W <x 2’ t)9 l (t)dt + c o I H y] (x 2' t)h l (t,dt + 7 I t^r g 2 (t)dt 


= -c Q (a y cos 2 9+a xx sin 2 e-a xy sin2e) , (c<x 2 <d) 


(40) 
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>21 


G x y( X 2’ t ) g i( t ) dt + C C 


H x] (x 2’ t)h l (t)dt + v 


t-x 


= -c o [(a~ y -a” x )sin9cose + a“ cos2e], (c<x 2 <d) 
l+K . _ >'o (k+1)(k o +,) 1 


v/nere 


c = 


'o 2y » G(x l } = 2u(< 0 -l) h^T 5 


H(x 1 ) = 


1 

2m h Q (x-j ) 


If there is no crack in the medium, g 2 =0=h 2 , the integral 
uncouple and (38) and (39) give the unknown functions g^ and h 
if the inclusion has an elliptic cross-section given by 


h 0 (x) = b Q /l-x z , 


(38) becomes 


1 


1 f f C 1 

V dt • j ^=? 9 i (t)dt = 

-i -i 1 x 


■c cr 

o yy 


where 


u 0 (1 + k)0 + k 0 ) 

C 1 2,b 0 ( v l) ’ 


and without any loss in generality it is assumed that a=-l, b= 
solution of (44) is found to be 

co 

9-1 (t) = -&• -±- • <-l<t<l) 

1 1 C 1 /PP 

which, for y Q =0 reduces to the well-known crack solution. By 
ing definition of the stress intensity factor 

k-j(t) = -lim ^2(l-x) g-j(x) , 

X-vl 


— h 9 (t)dt 
2 * 

( 41 ) 


(42) 

equations 
y For example, 

(43) 


(44) 

(45) 

1 , x^=x. The 

(46) 

using the follow- 

(47) 
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( 48 ) 


from (46) it follows that 

oo 

k m = 

V 1 ' 1+C-j * 


Similarly, in the absence of a crack from (39), (42) and (43) it may be 
shown that 


co 


h-i(t) = - — , (-l<t<l) 

1 1+c 2 /PF 

a xv v 0 (l +l <) 

ko(l) = lx. ’ C 9 = 


l+c 2 * 2 


2yb 


(49) 


(50) 


As another special case if we assume that the stiffness of the inclusion 
y Q =0, then the functions G and H defined by (42) vanish and the integral equa- 
tions ( 38) - (41 ) reduce to that of two arbitrarily oriented cracks shown in 
Fig. 1. 


3. Stress Intensity Factors 

In the linearly elastic medium under consideration the intensity of the 
stress state around the end points of the crack and the inclusion is governed 
by the singular behavior of the displacement derivatives g-j , g 2 , h-j and h 2 
which are defined by (l)-(4). If we assume the following standard definition 
of Modes I and II stress intensity factors 

k-| (a) = lim /2(a-x-j) CT yy( x -|>°) , (51) 

k 2 (a) = lim ^(a^) a] y (x.,,0) , (52) 

k.j(c) = lim /2(c-x 2 ) a yy (x 2 ,0) , etc. , (53) 

X2+C 

and observe that the system of integral equations ( 38)- ( 41 ) which has simple 
Cauchy type kernels has a solution of the form 
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( 54 ) 


G,(t) H - ( t ) 

g,(t) = 1 , h i (t) = 1 ■ ■■■ , (1=1.2) , 

/(b-t) (t-a) /(d-t)(t-c) 

from ( 38) -(41 ) and (51)-(54) it can be shown that 

k -j ( a ) = ^ lim /2(x.j-a) g-, (x-j ) , (55) 

X*| 

k.j(b) = - lini f2(b-x^ ) g-jU-j) , (56) 

x-j->b 

Ma) = 4=f“ lim /2(x-j-a) h^x-,) , (57) 

k«(b) = - lim /2(b-x^ ) h-^x-j) . (58) 

K x-,-*b 


The stress intensity factors (c) and k^(d), (i=l,2) may be expressed in 
terms of g£ and ^ ^ means of equations similar to ( 55)- ( 58) . 

4. Results 

The integral equations ( 38) - ( 41 ) are solved by using the technique 
described in [4] and the stress intensity factors are calculated from (55)- 
(58) and from similar expressions written for the crack. For various crack- 
inclusion geometries and stiffness ratios y Q /y (y Q being the shear modulus of 
the inclusion) the calculated results are given in Tables 1-6. The main 
interest in this paper is in relatively "thin" and flat inclusions. Hence 
in the numerical analysis it is assumed that the thickness h Q is constant. 
Table 1 shows the normalized stress intensity factors in a plane which con- 
tains a crack equal in size and coplanar with an inclusion and subjected to 
uniform tension and shear away from the crack-inclusion region (Fig. 2a). 

The inclusion model used in this analysis is basically a crack the surfaces 
of which are held together by an elastic medium of shear modulus y Q . Thus, 

for y =0 one recovers the two crack solution. It may be observed that for 
o 

y Q >0 there is a significant reduction in the stress intensity factors around 
the end points x-j=a and x^=b (Fig. 2a). In Table 1 the variables are the 
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stiffness ratio y /y and the thickness of the inclusion h Q /a-j with the 
spacing a/a-, = 0.01 being constant, where 2a-, is the length of the inclusion 
(Fiq. 2a). Similar results calculated by assuming that h Q /a^ = 1/20 and 

a/a-| is variable are shown in Table 2. 

For various values of the stiffness ratio y Q /y and fixed values of the 
inclusion thickness (h^a^l/20) and the distance a (a/a^O.l), the effect 
of the angle 0 on the crack tip stress intensity factors are given in Table 
3. The geometry and the loading condition away from the crack-inclusion 
region are shown in Fig. 2b. In this example, too, it is assumed that the 
inclusion and the crack are of equal length (a 2 =a-,). For the special case 
of y =0, that is, for the case of two cracks of equal lengths oriented at an 
angle e the stress intensity factors are given in Table 4. 

The stress intensity factors for the symmetric crack-inclusion geometries 
shown in Figures 3a and 3b are given in Table 5, where the length ratio a 2 /a-, 
is assumed to be the variable. In both examples the inclusion (half) length 
a is used as the normalizing length parameter and the relative distance 
c/a-, (Fig. 3a) or a/a-j (Fig. 3b) is assumed to be constant. 

Table 6 gives the stress intensity factors for a crack perpendicular 
to the inclusion where, . referring to Fig. 1, 9=tt/2, a=0, y Q -y/20 and c/a-,- 
0.05* are fixed and a 2 is variable. 

It should be noted that since the superposition is valid, the tables 
give the stress intensity factors for the most general homogeneous loading 
conditions away from the crack-inclusion region. Also, the tables give the 
stress intensity factors which are normalized with respect to where 2a-, 

is the length of the inclusion and (i,j)=(x,y), (Fig. 1). The notation used 
in the tables is 


_ "1 


k,(a) 


'la 




, k 


2a 


k 2 (a) 


, k 


lc 


_ 1 


Me) 




, etc. 


(59) 


where k-, and k 2 are, respectively. Modes I and II stress intensity factors 
defined by equations such as (51 )-(53) and calculated from the expressions 

such as (55)-(58). 


- 11 - 



REFERENCES 


1. F. Erdogan and G.D. Gupta, "The inclusion problem with a crack crossing 
the boundary", Int. Journal of Fracture, Vol. 11, pp. 13-27, 1975. 

2. F. Erdogan, G.D. Gupta and M. Ratwani , "Interaction between a circular 
inclusion and an arbitrarily oriented crack", Journal of Applied Mechanics, 
Vol. 41, Trans. ASME, pp. 1007-1013, 1974. 

3. J. Dundurs, "Elastic interaction of dislocations with inhomogeneities". 
Mathematical Theory of Dislocations , T. Mura, ed., pp. 70-115, ASME, 

New York, 1969. 

4. F. Erdogan, "Mixed boundary value problems in Mechanics", Mechanics 
Today, S. Nemat-Nasser, ed. , Vol. 4, pp. 1-81 , 1978. 


- 12 - 



Table 1. Modes I and II stress intensity factors for the case of a crack 
located in the plane of the inclusion in a medium subjected to 
a°°, or cTy away from the crack-inclusion region (Fig. 2a); c=-a, 
d=-b, a/a-pO.Ol , k-j c =ki (c)/Cyy' / aY, k-j^k-j (d)/ayyV / a!jr, k 2 C ~k 2 (c)/cr X yV / iY’ 

k 2d =k 2^ d ^ /,CT xy'^T 5 k la = M a ^ a yy' /a 7’ k 2a =k 2^ a ^ a xy’ /a T’ k lb~M b ^ a yy v/a 7’ 
k 2b =k 2^ b ^ a xy y/a T * a i = ( b_a )/ 2 - 


2h o 



b-a 

0 


071 '” 

isai 

075 

- T70 

O 

570 — 


0.01 

1.2063 

.1578 

.1031 

.0535 

.0303 

.0163 

.0085 


k.. 

0.02 

1.2063 

.2320 

.1578 

.0888 

.0535 

.0303 

.0163 


lb 

0.1 

1.2063 

.5146 

.3713 

.2320 

.1578 

.1031 

.0634 

.0303 


0.2 

1.2063 

.6836 

.5146 

.3323 

.2320 

.1578 

.1031 

.0535 


0.01 

2.9642 

.5725 


.2104 

.1207 

.0654 

■BMiKM 

.0140 

k. 

0.02 

2.9642 

.7941 


.3404 

.2104 

.1207 

.0654 

.0276 

la 

0.1 

2.9642 

1.5036 

1.1620 

.7941 

.5725 

.3908 

.2478 

.1207 


0.2 

2.9642 

1.8803 

1.5036 

1.0636 

.7941 

.5725 


.2104 


0.01 

2.9642 

1.1795 

1.1045 

1.0479 


1.0132 

1.0067 

1.0027 

k. 

0.02 

2.9642 

1.2952 

1.1795 

1.0870 

1.0480 

1.0255 

1.0132 

1.0054 

1c 

0.1 

2.9642 

1 . 7825 

1.5321 

1.2952 

1.1795 

1.1045 

1.0583 

1.0255 


0.2 

2.9642 

2.0764 

1.7825 

1.4645 

1.2952 

1.1795 

1.1045 

1.0479 


0.01 

1.2063 

1.0116 

1.0063 

HBtjsfcj 

1.0014 

1.0007 

1.0004 

1.0001 

i. 

0.02 

1.2063 

1.0211 

1.0116 


1.0027 

1.0014 

1.0007 

1.0003 

k ld 

0.1 

1.2063 

1.0693 

1.0432 


1.0116 

1.0063 • 

1.0033 

1.0014 


0.2 

1.2063 

1.1019 

1.0693 

in 

1.0211 

1.0116 

1.0063 

1.0027 


0.01 

1.2063 


.2159 

■1KH 

.0810 1 

mwm 

.0269 

.0117 

k OK 

0.02 

1.2063 


JmSEm 

ll I 

.1275 

■ 

.0482 


2b 

0.1 

1 . 2063 


SI 

.4368 

.3106 

.2159 

.1459 

.0810 


0.2 

1.2063 

.9673 

.8214 

.5946 

.4368 

.3106 

.2159 

.1275 


TJ70T - 

VSKIiGHI 

umm 


.4743 

mam u 

wmm 

.1076 

■ SfiuS 

k. 

0.02 

2.9642 

1.3214 


.6747 

.4743 

wsBm 



2a 

0.1 

2.9642 

2.1749 

iHkimI 

1.3214 

1.0075 

.7480 

.5345 

.3122 


0.2 

2.9642 

2.4785 

2.1749 

1.6847 

1.3214 


.7480 

.4743 


0.01 

2.9642 

1.4272 


1.1366 

1.0778 

1.0425 

1.0225 

1.0093 

k. 

0.02 

2.9642 

1.6463 

EK wm 

1.2298 

1.1366 

1.0778 

1.0425 

1.0182 

2c 

0.1 

2.9642 

2.3136 

2.0183 

1.6463 

1.4272 

1.2691 

1.1622 

1.0778 


0.2 

2.9642 

2.5619 

2.3136 

1.9221 

1.6463 

1.4272 

1.2691 

1.1366 


0.01 

1.2063 

1.0330 

1.0188 

1 . 0085 

1.0045 

1.0023 

1.0012 

1.0005 

k^ » 

0.02 

1.2063 

1.0549 

1.0330 

1.0156 

1.0085 

1.0045 

1.0023 

1.0010 

2d 

0.1 

1.2063 

1.1292 

1.0954 

1.0549 

1.0330 

1.0188 

1.0103 

1.0045 


0.2 

1.2063 

1.1583 

1.1292 

1.0846 

1.0549 

1.0330 

1.0188 

1.0085 
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Table 2. Modes I and II stress intensity factors for the case of a crack 
located in the plane of the inclusion in a medium subjected to 
a yy ° r a xy away from the crack ‘ 1nc ' lus1 ' on re 9ion (Fig. 2a); c=-a, 
d=-b, h o /a 1 ®l/20. 



2a 

P 

O 

p. 


b-a 

0 

0.0b 

mam 

0.25 

0.5 i 

— 

2.0 

5.0 


0.01 

1.2063 

.3713 

.2611 

.1578 

.1031 

.0635 

.0366 

.0163 

k i . 

0.5 

1.0517 

.3544 

.2513 

.1527 

.0998 

.0615 

.0354 

.0158 

lb 

1 

1.0280 

.3493 

.2479 

.1508 

.0986 

.0607 

.0350 

.0156 


2 

1.0125 

.3453 

.2452 

.1492 

.0976 

.0601 

.0347 

.0154 


0.01 

2.9642 

1.1620 

.8751 

.5725 


.2478 

.1454 

.0654 

k. 

0.5 

1.1125 

.3877 

.2768 

.1693 

.1110 

.0685 

.0395 

.0176 

1 a 

1 

1.0480 

.3604 

.2564 

.1563 

.1023 

.0630 

.0364 

.0162 


2 

1.0176 

.3481 

.2474 

.1506 

.0985 

.0607 

.0350 

.0156 

■ 

0.01 

2.9642 

1.5321 

1.3433 

1.1795 

1.1045 

1.0583 

1.0313 


k. 

0.5 

1.1125 

1.0229 

1.0130 

1.0057 

1.0030 

1.0015 

1.0008 

1.0003 

1c 

1 

1.0480 

1.0096 

1.0054 

1.0024 

1.0012 

1.0006 

1.0003 

1.0001 


2 

1.0176 

1.0035 

1.0020 

1.0009 

1.0004 

1.0002 

1.0001 

1.0000 


0.01 

1.2063 

1.0432 

1.0253 

1.0116 



IBSiSUI 

■Hi i iflM 

k, . 

0.5 

1.0517 

1.0104 

1.0058 

1.0026 

1 


HS9 


id 

1 

1.0280 

1.0056 

1.0031 

1.0014 



1.0002 

1.0001 


2 

1.0125 

1.0025 

1.0014 

1.0006 



1.0001 

1.0000 



1.2063 

.6500 

■Hn 

.3106 

.2159 

.1459 

.0943 


k n . 

0.5 

1.0517 

.6031 

mWlm 

.2979 


.1412 

.0914 

.0467 

2b 

1 

1.0280 

.5925 

.4503 

.2938 

.2057 

.1395 

.0903 

.0461 


2 

1.0125 

.5849 


.2905 

.2035 


.0893 

.0456 


0.01 

2.9642 , 

1.8071 

1.4340 

1.0075 j 


.5345 

.3601 


k« a 

0.5 

1.1125 

.6498 

.4971 

.3272 , 

.2302 

.1567 

.1017 


2a 

1 

1.0480 

.6081 

.4636 

.3035 j 

.2129 

.1446 

.0937 



2 

1.0176 

.5889 

.4483 

.2930 

.2053 

.1393 




0.01 

2.9642 

2.0183 

1.7299 

1.4272 

1.2691 

1.1623 

j 1 .0937 

1.0425 

k o _ 

0.5 

1.1125 

1.0523 

1.0344 

1.0172 

1.0095 




2c 

1 

1.0480 

1.0222 

1.0145 

1.0072 

1.0040 





2 

1.0176 

1.0081 

1.0053 

1.0026 

1.0014 





0.01 

1.2063 

1.0954 

1.0637 

1 .0330 

1.0188 

1.0104 

1.0055 

1.0023 


0.5 

1.0517 

1.0239 

1.0157 

1.0078 

1.0043 

1.0023 

1.0012 

1.0005 

k 2d 

1 

1.0280 

1.0129 

1.0084 

1.0042 

1.0023 

1.0012 

1.0006 

1.0003 


2 

1.0125 

1.0057 

1.0038 

1.0019 

1.0010 

1.0005 

1.0003 

1.0001 
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Table 3. 


The effect of angular orientation a and the modulus ratio 
uj/y on the stress intensity factors in a medium under general 
in-plane loading (Fig. 1); c=a, d=b, 2h o /(b-a)=l/20, 2a/(b-a)=0.1 . 


i 1 

1, 

9 w 

I 3 j 

K 

30 

60 

90 | 

120 1 

150 i 

180 

r 




y o /y = 

0.05 




k lc 

0.2624 

0.8047 

1.0961 

0.8097 

0.2654 

0 

00 

CT 

k 2c 

-0.4711 

-0.4636 

0.0163 

0.4737 

0.4585 

0 

XX 

k ld 

0.2560 

0.7618 

1.0106 

0.7562 

0.2518 

0 


k 2d 

-0.4378 

-0.4253 

0.0122 

0.4432 

0.4383 

0 


C 

0.6402 

0.2232 

-0.0311 

0.2749 

0.8366 

1.1094 

CO 

k 2c 

0.4596 

0.4217 

-0.0483 

-0.5019 

-0.4771 

0 

c yy 

k ld 

0.7052 

0.2221 

-0.0109 

0.2568 

0.7702 

1.0250 


k 2d 

0.4105 

0.3981 

-0.0386 

-0.4636 

-0.4493 

0 


k lc 

-0.5020 

-0.5895 

0.2839 

1.1440 

1.0302 

0 

CO 

CT , 

k 2c 

0.3394 

-0.5681 

-1.0010 

-0.3793 

0.7098 

1.2367 

xy 

k ld 

-0.9072 

-0.8566 

0.0354 

0.9049 

0.8903 

0 


k 2d 

0.4353 

-0.5284 

-0.9911 

-0.4631 

0.5521 

1.0567 

1 


— r 1 

V y = 

0.1 

♦ 


k lc 

0.2552 

0.7786 

1.0613 

0.7908 

0.2608 

0 

oo 

k 2c 

-0.4593 

-0.4546 

0.0095 

0.4610 

0.4512 

0 

a xx 

k l d 

0.2534 

0.7570 

1.0066 

0.7540 

0.2512 

o 1 

i 


k 2d 

-0.4366 

-0.4291 

0.0072 

0.4395 

0.4366 

0 


k 2c 

0.6535 

0.2334 

-0.0181 

0.2628 

0.8003 

1.0643 


k 0 

0.4533 

0.4238 

-0.0293 

-0.4758 

-0.4605 

0 

CO 

2c 







a yy 

k ld 

0.7248 

0.2350 

-0.0058 

0.2540 

0.7615 

1.0143 


k 2d 

0.4219 

0.4145 

-0.0215 

-0.4506 

-0.4425 

0 


k lc 

-0.6023 

-0.6717 

0.1849 

1 . 0482 

0.9749 

0 

oo 

k 0 

0.3956 

-0.5401 

-0.9996 

-0.4197 

0.6414 

1.1599 

0„ w 

2c 







xy 

k l d 

-0.8892 

-0.8588 

0.0230 

0.8910 

0.8817 

0 


k 2d 

0.4617 

-0.5172 

-0.9943 

-0.4762 

0.5343 

1.0374 
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Table 3 - cont. 





V y 

= 0.5 




k lc 

0.2478 

0.7537 

1.0157 

0. 7622 

0.2535 

0 

CO 

a vv 

k 2c 

-0.4414 

-0.4405 

0.0019 

0.4418 

0.4391 

0 

XX 

k ld 

0.2509 

0.7517 

1.0017 

0.7511 

0.2503 

0 


k 2d 

-0.4341 

-0.4322 

0.0017 

0.4347 

0.4341 

0 


k lc 

0.7013 

0.2427 

-0.0045 

0.2523 

0.7620 

1.0158 

00 

k 2c 

0.4381 

0.4288 

-0.0078 

-0.4448 

-0.4407 

0 

a yy 

k ld 

0.7446 

0.2469 

-0.0012 

0.2510 

0.7527 

1.0033 


k 2d 

0.4312 

0.4292 

-0.0048 

-0.4371 

-0.4353 

0 


k lC 

-0.7657 

-0.8011 

0.0517 

0.9166 

0.3971 

0 


k 2c 

0.4738 

-0.5057 

-0.9981 

-0.4766 

0.5420 

1.0479 


k ld 

-0.8712 

-0.8639 

0.0061 

0.8726 

0.8702 

0 


k 2d 

0.4910 

-0.5046 

-0.9987 

-0.4938 

0.5094 

1.0105 




• 

y o /y 

= 2 




k lc 

0.2484 

0.7504 

1.0041 

0.7535 

0.2510 

0 

CO 

a vv 

k 2c 

-0.4356 

-0.4354 

0.0003 

0.4356 

0.4349 1 

0 

XX 

k ld 

0.2503 

0.7505 

1.0004 

0.7503 

0.2501 

0 


k 2d 

-0.4333 

-0.4328 

0.0004 

0.4335 

0.4333 

0 


k lc 

0.7317 

0.2473 

-0.0012 

0.2505 

0.7531 

1.0042 

00 

k 2c 

0.4330 

0.4314 

-0.0022 

-0.4363 

-0.4352 

0 

yy 

k ld 

0.7487 

0.2492 

-0.0003 

0.2503 j 

0.7507 

1.0009 


k 2d 

0.4326 

0.4321 

-0.0012 

*3- 

CO 

• 

o 

1 

-0.4336 

0 


k lc 

-0.8318 

-0.8460 

0.0146 

0.8801 

0.8748 

0 

CO 

a 

k 2c 

0.4947 

-0.5001 

-0.9989 

-0.4932 

0.5122 

1.0139 

xy 

k ld 

-0o8674 

-0.8655 

0.0016 

0.8678 

0.8672 

0 


k 2d 

0.4976 

-0.5013 

-0.9997 

-0.4983 

0.5026 

1.0029 
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Table 4. Interaction of two cracks (Fig. 2b); Vq/u-j® 0 , c-a, d-b, 
2a/(b-a) = 0.1. 



e° 

30 

60 

90 

120 

150 

180 


k la 

0.1834 

-0.0122 

-0.1604 

-0.1271 

-0.0361 

0 


k 2a 

0.1293 

0.0928 

0.2122 

0.2877 

0.1946 

0 


k ib 

-0.1471 

-0.1373 

-0.0566 

-0.0113 

0.0024 

0 

00 

k 2b 

0.1825 

0.2323 

0.2104 

0.1371 

0.0588 

0 

a xx 

k lc 

0.3637 

1.0032 

1.2370 

0.8684 

0.2790 

0 


k 2c 

-0.5576 

-0.4950 

0.0577 

0.5191 

0.4810 

0 


k id 

0.3073 

0.8057 

1.0308 

0.7633 

0.2536 

0 


k 2d 

-0.3956 

-0.3708 

0.0477 

0.4591 

0.4441 

mm 


^la 

0.5843 

0.9140 

1.2370 

1.3954 

1.4643 

1 .4914 


k 2a 

-0.1912 

-0.0242 

-0.0577 

-0.1080 


0 


k lb 

0.9210 

1.0081 

1.0308 


1.0994 

1.1220 

oo 

k 2b 

0.0215 

-0.0427 

-0.0477 

-0.0168 


0 

a 

yy 

k lc 

0.4051 

-0.1004 

-0.1604 

0.3999 

1.1491 

1.4914 


^2c 

0.6195 

0.4264 

| 

-0.2122 

-0.6987 

-0.6027 

0 


k id 

0.4666 

0.0652 

-0.0666 

0.2821 

0.8481 

1 .1220 


k 2d 

0.1916 

0.1811 

-0.2104 

-0.5795 

-0.5082 

0 


m 


0.7402 

0.6381 

0.3381 

0.1384 

0 


IS 


1.1315 

1.0152 

1.1777 

1.4058 

1.4914 


k lb 

0.4327 

0.1938 

0.0748 

0.0610 

0.0532 

0 

oo 

k 2b 

0.5851 

0.7960 

0.9950 

1.1104 

1.1305 

1.1220 

a xy 

k lc 

-0.4402 

-0.4311 

0.6381 

1.4876 

1.2302 

0 


k 2c 

0.3095 

-0.6671 

-1.0152 

-0.2462 

0.9347 

1.4914 


k ld 

-1.1414 

-0.8951 

0.0748 

0.9554 

0.9234 

0 


k 2d 


-0.6362 

-0.9950 

-0.4219 

0.6115 

1.1220 
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Table 5. Stress intensity factors for the case of a crack perpendicular 
to the inclusion, y o /u=l/20, h o /a-|=l/20. 


Fig. 3a 
a=-b=-a-j 
c/a-|=0. 1 

CO 

rr 

k 

a 2 /a l 

0.1 

0.5 

1.0 

5.0 

CO 

a xx 

k,,-k n 

1 U |Lf 

k 2a = " k 2b 

k lc 

k ld 

-0.0088 

-0.0058 

+1.0536 

1.0320 

-0.0479 

-0.0820 

1.1611 

1.0245 

-0.0938 

-0.1428 

1.1572 

1.0109 

-0.1449 

-0.2729 

1.1256 

1.0029 

CO 

0 

yy 

k la =k lb 

k 2a =_k 2b 

k lc 

k ld 

0.3424 

0.0006 

-0.1220 

-0.0988 

0.3441 

0.0039 

-0.0896 

-0.0116 

0.3441 

0.0039 

-0.0632 

0.0067 

0.3438 

0.0033 

-0.0255 

0.0021 

CO 

a xy 

k la = " k lb 

k 2a =k 2b 

k 2c 

k 2d 

-0.0004 

0.5703 

-0.7288 

-0.7856 

-0.0162 

0.5162 

-0.9533 

-1.0338 

-0.0850 
0.4502 
-1 .0730 
-1.0638 

-0.5164 

0.4199 

-1.2431 

-1.0200 

Fig. 3b 
c=-d=-a2 
a/a-j =0. 1 

CO 

a xx 

k la 

k lb 

k lc =k ld 

k 2c =_k 2d 

0.0208 

0.0006 

1.0037 

-0.0011 

-0.1238 

0.0100 

1.0053 

-0.0074 

-0.2149 

0.0234 

1.0101 

-0.0107 

-0.2773 

-0.1170 

1.0026 

-0.0045 

CO 

a yy 

k la 

k lb 

k lc =k ld 

k 2c =_k 2d 

0,3476 

0.3416 

0.1584 

-0.0353 

0.3543 

0.3418 

-0.0186 

0.0460 

0.3764 

0.3416 

-0.0324 

0.0406 


00 

a xy 

k 2a 

k 2b 

k lc = - k ld 

k 2c =k 2d 

0.6514 

0.5808 

-0.4813 

-1.3694 

0.5903 

0.6066 

-0.2431 

-0.9632 

0.4304 
C. 631 5 
-0.1012 
-0.9372 

0.0544 

0.3702 

-0.0010 

-0.9946 
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Table 6. Stress intensity factors for a crack perpendicular to the 
inclusion (Fig. 1); 9 =t/ 2, a=0, 2c/(b-a)=Ch05, u /u=l/20, 
2h Q /(b-a)=0.05. 


00 

b 


V a i 



a 

K 

0.1 

0.5 

1.0 

5.0 


^la 

.0399 

.2055 

.3675 

1.1277 


^2a 

.0128 

.0418 

.0555 

.1125 


k lb 




-.0715 

oo 

k 2b 


.0402 


.3050 

a xx 

k lc 

1.0762 

1.1674 

1.1729 

1.1435 


^2c 

.0162 

-.0056 




k ld 

1.0310 

1.0274 

1.0143 



k 2d 

.0207 

.0212 

.0115 



k la 

7sm 

737T5 


.3884 


k 2a 

.0092 

.0283 

.0390 

.0533 


k lb 

.3414 

.3411 

.3418 

.3456 

00 

k 2b 

.0001 

-.0010 

-.0036 

-.0062 

a yy 

k lc 

-.0490 

-.0607 

-.0514 

-.0250 


k 2c 

-.3157 

-.2298 

-.1863 

-.0933 


k ld 

-.0468 

-.0250 

-.0084 

-.0009 


k 2d 

-.1943 

-.0830 

-.0464 

-.0048 


m 

.0887 

.3231 

.4952 

1.1795 



.6265 

.7947 

.9710 

1.9112 


— 

.0002 

.0001 

.0079 

.2709 

CO 

■9 

.5805 

.5910 

.5713 

.4743 

a xy 

k lc 

1.1620 

.6411 

.4373 

.1825 


k 2c 

-1.0423 

-1.1380 

-1.1889 

-1.2670 


k ld 

.6504 

.1454 

.0426 

.0045 


k 2d 

-.9292 

-.9710 

-1.0075 

-1.0117 
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